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ABSTRACT 

An abundance decrease in carbon- and oxygen-bearing species relative to dust has been fre¬ 
quently found in planet-forming disks, which can be attributed to an overall reduction of gas 
mass. However, in the case of TW Hya, the only disk with gas mass measured directly with 
HD rotational lines, the inferred gas mass (<0.005 solar mass) is significantly below the directly 
measured value (^0.05 solar mass). We show that this apparent conflict can be resolved if the 
elemental abundances of carbon and oxygen are reduced in the upper layers of the outer disk but 
are normal elsewhere (except for a possible enhancement of their abundances in the inner disk). 
The implication is that in the outer disk, the main reservoir of the volatiles (CO, water, ...) 
resides close to the midplane, locked up inside solid bodies that are too heavy to be transported 
back to the atmosphere by turbulence. An enhancement in the carbon and oxygen abundances 
in the inner disk can be caused by inward migration of these solid bodies. This is consistent with 
estimates based on previous models of dust grain dynamics. Indirect measurements of the disk 
gas mass and disk structure from species such as CO will thus be intertwined with the evolution 
of dust grains, and possibly also with the formation of planetesimals. 

Subject headings: astrochemistry — circumstellar matter — molecular processes — planetary systems 
— planet-disk interactions — planets and satellites: atmospheres 


1. Introduction 


Observations and modeling of transitional disks 
have pointed to an apparent lack of gas relative 
to their dust content, particularly at larger dis¬ 
tances from the star. Early observations (see, 
e.g. [Skrutskie et al.||1991| and [Yamashita et al?] 
1993) show that the CO emission are much lower 


than expected from dust emission and molecular 
cloud composition. A recent Herschel/PACS sur¬ 
vey finds that transitional disks have weaker O l 
63 jam emission than their full disk counterparts, 
even if they have similar continuum at 63 jam 


(Keane et al. 2014 Howard et al. 2013). One 


interpretation of this is that the gas mass is re¬ 
duced, or, equivalently in effect, that the transi¬ 
tional disks are less flared than their full disk coun¬ 


terpart (Keane et al. 2014). Aresu et al. (2014) 


find that with a normal gas-to-dust ratio, mod¬ 


els of the O I 63 jam emission for disks in Taurus 
tend to over-predict the line intensity. Similarly, 
detailed modeling of the prototypical transitional 
disk, TW Hya, suggest a small gas mass (0.5- 
5)xl0“^ Mq (and small gas-to-dust mass ratio) 
(Thi et al. |2010[ Kamp et al. 2013 
Best||2014 ). 


Williams & 


However, the gas mass in these studies are 
only inferred rather than directly measured. The 
main constituent of the gas is hydrogen, mostly in 
molecular form, followed by helium, both of which 
defy direct detection in the bulk of the disk, be¬ 
cause their line excitation requires high tempera¬ 
ture. So the gas mass is inferred from other agents, 
such as the commonly used CO lines or dust con¬ 
tinuum, in combination with radiative transfer 
and thermo-chemical models. Fortunately, with 
Herschel, it is possible to measure the gas mass 
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in a much more direct way, by observing the ro¬ 
tational transition lines of HD, an isotopologue of 
molecular hydrogen. The translation of HD inten¬ 
sity to H 2 mass is more straightforward because 
there is no complex chemistry involved (though 
the result depends on the temperature structure). 
Direct measurement of the gas mass of TW Hya 
based on Herschel observation of HD lines gives 
mass ^0.05 Mq ( Bergin et al.|2Q13 ), much higher 
than the values in some of the previous studies. 

To reconcile the two classes of seemingly dis¬ 
parate results on the gas mass of TW Hya, here 
we show that it is not the overall gas mass — the 
mass of atomic and molecular hydrogen — that 
is reduced, but rather it is the mass of the ice¬ 
forming species, such as CO and H 2 O, are reduced 
in the upper layers of the outer disk. It is these 
upper layers that are emissive and detectable by 
astronomical observations. The depletion of CO 
relative to hydrogen in DM Tan and GG Tan has 


been found by Dutrey et al. (1997) based on exci¬ 


tation arguments, and in TW Hya by |Favre et al. 
(2013) based on direct comparison of the CO iso¬ 
topologue emission with the HD result. Reduc¬ 
tion of the water content has been suggested by 
Bergin et al. I ( |2Q1Q| ) and [Hogerheijde et ah] ( |2Q11[ 
supporting material). Since CO and H 2 O are the 
main bearers of carbon and oxygen in the upper 
layers (below the photodissociated layer), the to¬ 
tal abundances of elemental carbon and oxygen 
are lowered. This affect other species chemically 
linked to them. We also establish that this effect 
is evident only in the outer disk (beyond the water 
and CO snow line), while in the inner disk the op¬ 
posite — namely enhancement of the oxygen and 
carbon abundances — might be taking place. A 
possible physical mechanism based on dust grain 
settling and migration is provided at the end. 

2. Models 

We run two models for TW Hya using a thermo¬ 


chemical and radiative-transfer code (Du & Bergin 
2014). The system is evolved for 1 Myr. Both 


models assume a disk with total gas mass 0.05 Mq 
as constrained by the HD rotational lines (Bergin 
a distance of 51 pc ( |Mamajek]|2005[) , 
and an inclination angle of 7 degree (Qi et al. 


2004). The dust mass is not very well constrained. 


for both models. We assume three components for 
the disk structure: an optically thin component in¬ 
side ^3.5 AU, a component tapering at ^50 AU 


(Andrews et al. 2012 


Menu et al. 


2014), and a 


component extended to ^200 AU. The mass ra¬ 
tio of the three components are 10“^ : 20 : 1. 
The outer two components have a surface density 
profile S oc The disk 3-D density and tem¬ 

perature structure are determined based on ver¬ 
tical hydrostatic equilibrium, calculated through 
iterative Monte Carlo simulation of dust radiative 
transfer, which also yields a spectral energy dis¬ 
tribution that matches observation. The result¬ 
ing 870 /im continuum map of TW Hya resem¬ 
bles the observed one of Andrews et al. (2012), 


though we spend no effort to match the visibility 
profile. The distribution of gas density and tem¬ 
perature is shown in panel (a) and (b) of Fig. 
The only difference between the two models lies 
in the abundances of oxygen and carbon (relative 
to hydrogen). The first model has canonical “low- 
metal” ISM (interstellar medium) oxygen and car¬ 
bon abundances ( [Garrod et al.|2008 ) in the whole 
disk, while the second model has reduced oxygen 
and carbon abundances in the upper layers of the 
outer disk, and enhanced oxygen and carbon abun¬ 
dances in the inner disk (see panel (c) and (d) of 
Fig. [^. The degree of reduction for oxygen and 
carbon are based on an analytical prescription, the 
exact form of which is manually adjusted for a best 
fit. Carbon is less reduced than oxygen in the 
outer disk, inspired b y the difference in the evap¬ 
oration temperatures ( Oberg et al.||2Qll ) of water 
and the major bearers of carbon (CO and CO 2 ). 
The enhancement in the inner disk is simply im¬ 
plemented with a step function. This empirical 
approach is particularly tailored to match the ob¬ 
servational data. Future self-consistent study is 
needed. 

3. Results 

The code predict intensities for a suite of transi¬ 
tions of different species, including C ll, O l, H 2 O, 
HD, and the CO isotopologues, which have been 
detected by observation. The resulting line inten¬ 
sities from the two models are compared with the 
observed values in Fig. (for the numerical values 
see Table . In agreement with previous stud- 


so we use a uniform gas-to-dust mass ratio of 100 


ies (Thi et al. 2010[ |Kamp et ah 2013 Williams 


fc Best||2014 Keane et al.||2014' ), the model with- 
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Fig. 1.— Distribution of disk physical parameters 
used in the models, (a): gas density distribution 
in the two models, (b): gas temperature distri¬ 
bution in the model in which the oxygen and car¬ 
bon abundances are changed; for the model with 
full oxygen and carbon abundances, the gas tem¬ 
perature distribution is slightly different due to 
changes in the heating and cooling rates of oxygen- 
and carbon- bearing species, (c) and (d): distribu¬ 
tion of oxygen and carbon abundances relative to 
hydrogen nuclei in the model in which their abun¬ 
dances are changed. 


out elemental depletion consistently over-predicts 
the line intensities of CO, water, and a few other 
species for TW Hya, while the model with depleted 
oxygen and carbon abundances matches the ob¬ 
servational data. There are a few lines that most 
clearly require oxygen and carbon to be depleted 
in the atmosphere of the outer disk, including the 
O I 63 /im, the CO and C^^O low-J lines, and all 
the water ground state lines. 
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Fig. 3.— Fitting to the Spitzer/IRS H 2 O and 
OH lines. The Spitzer spectrum (in gray) is from 
Zhang et al. (2013). 


In both models, the water lio — loi? lii ~ Ooo? 
and 2 12 — loi lines originate from >20 AU, well 
beyond the water snow line. Gas phase wa¬ 
ter molecules in this region are mainly produced 


through photodesorption of water ice (Bergin et al. 

2010 

Hogerheijde et al.||2011 Kamp et al. 

20131), 

with their rates established by experiment ( 

Oberg 


et al. 2009). A normal oxygen abundance in the 


outer disk over-predicts their intensities by more 
than one order of magnitude. The situation is 
similar for the CO and C^^O low-J lines (see 

The O I and C ll lines 


Favre et al. 


2013) 


also 

mainly originate from near the water snow line 
at r < 4 AU in the model with reduced oxygen 
and carbon abundances, coincident with where the 
water infrared lines originate ( Zhang et al.||2013 ). 
With normal oxygen and carbon abundances, ad¬ 
ditional contribution to the emission arises in gas 
out to 10 - 20 AU, which make them too strong, 
by a factor of 5 in the case of O l. We note that 
a deficit of elemental carbon in the disk atmo¬ 
sphere has been suggested for the Herbig Be star 
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Fig. 2.— Ratios between the modeled and observed intensities for a selection of lines. Red dot: with full 
oxygen and carbon abundances; blue rectangle: with depleted oxygen and carbon abundances. The arrows 
mean that the observed values are upper limits. The relative measurement error of the observational data 
is usually much less than 50%. For the detected lines, a perfect fit would land on the dashed line. 
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Table 1: Comparison between modeled and observed integrated line fluxes. The “Full” column is for the 
model with canonical carbon and oxygen abundances, and the “Depleted” column is for the model with 
reduced carbon and oxygen abundances. 


References. — 1. ( [Bergin et al.||2013|>; 2. (|Favre et al.||2013|>; 3. ([Hogerheijde et al.||2011|); 4. jQi et al.||2006| >; synthesized 
beam corresponding to 100 AU 5. ( |Thi et al.|2010| ); 6. ( [Zhang et al.|2013| >. 7. ( [Tsnkagoshi et al.|2015| >. 
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HD 100546 by Bmderer et al. (2012). 

In the model with elemental depletion, we let 
the carbon and oxygen abundances gradually re¬ 
turn to normal when approaching the CO snow 
line (^20 AU) from outside. This is required to 
account for the observed high-J and rovibrational 
CO lines, as well as O l and the numerous mid- 
infrared transitions of water vapor (see Fig. [^. 
This scenario resonates with the recent result on 
the CO/H 2 ratio, which is ^ 1.6 x 10 “^ very close 
to the star in RW Aurigae A ( France et al.||2Q14 ). 
For the ^^CO 10—9 and CO 10—9 and 23—22 tran¬ 
sitions, canonical carbon and oxygen abundances 
in the inner disk of the depleted model are not 
enough to reproduce the data; these three lines are 
under-predicted by about one order of magnitude. 
Along our line of reasoning, one intriguing possi¬ 
bility is that the carbon and oxygen abundances 
are enhanced (van Dishoeck et al. 2014), rather 
than reduced, in the inner disk. In fact, if we in¬ 
crease the abundance of carbon and oxygen by a 
factor of 10 over the ISM value within ^15 AU (see 
the insets in panel (c) and (d) of Fig. [^, the CO 
10 — 9 and 23 — 22 and the ^^CO 10 — 9 lines would 
be in good agreement with observations, as shown 
in Fig. while still keeping other lines consistent 
with observations. Inward migration and evapo¬ 
ration of icy bodies is predicted to cause such an 
enhancement of CO and water vapor in the inner 


disk ( Cnzzi fc Zahnle||2004 Ciesla fc Cnzzi||2006 ) 
In the above we have assumed that only oxy¬ 
gen and carbon are differentially depleted in the 
outer disk atmosphere, and have left nitrogen in¬ 
tact. This is a reasonable assumption, considering 
that solar system ices as traced by comets are well 
known to be nitrogen-poor ( Wyckoff et al.|[T99I ). 
One interesting outcome of this assumption is that 
by depleting oxygen (and carbon, but to a lower 
degree) over nearly the entire disk the gas will be¬ 
come rich in carbon (particularly inside the CO 
snow line) and nitrogen-bearing molecules. As an 
example, Fig.|^shows that the column densities of 
species containing C-N, C-H, and N-H bonds can 
be increased by up to two orders of magnitude in 
the depleted model. This is in line with the re¬ 
cent results on the enhanced abundances of C 2 H 


and CN in evolved protoplanetary disks (Kastner 


Depletion of carbon and oxygen in the outer 
disk also affect the HD emission. This is because 


the cooling efficiency due to species bearing these 
two elements are changed, so that the tempera¬ 
ture in the outer disk becomes higher, leading to 
a stronger HD emission. For example, 80% of the 
HD 112 jam emission comes from within 80 AU in 
the model with no depletion, while with depletion 
60% arises within the same radius (i.e. the outer 
disk contributes more). Though, as pointed out 
by Rollig et al. ( 2007[ ), the energy balance calcula¬ 
tion could be subject to considerable uncertainty 
in the tenuous regions. 



Fig. 4.— Comparison between the column den¬ 
sity of non-oxygen-bearing cyanide (containing C- 
N bonds), hydrocarbon (containing C-H bonds), 
and hydronitrogen (containing N-H bonds) from 
two models with full (solid lines) and reduced 
(dashed lines) oxygen and carbon abundances. 


4. Discussions 


The depletion of oxygen and carbon relative 
to hydrogen solves the puzzle regarding the gas 
mass of the TW Hya disk. The only known gas- 
loss mechanism that preferentially removes species 
containing oxygen and/or carbon (mainly H 2 O 
and CO) rather than atomic or molecular hydro¬ 
gen is ice mantle formation on dust grains. We 
emphasize that ice mantle formation and grain 
surface chemical reactions by themselves do not 
reduce the local elemental abundances (gas + ice) 
in the atmosphere. Reboussin et al.| ( |2Q15[ ) find 


that CO can be transformed into other species 
such as CO 2 and CH 4 on the dust grain surfaces 
(see also Bergin et al. 2014). However, what we 
find is that if the ice-coated grains stay in the 
upper atmosphere, desorption and dissociation by 
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UV photons will release H 2 O, O, OH, CO, etc. to 


years of the condensation of calcium-aluminum- 


the gas phase, and the emission from them will be 

rich inclusions (CAIs) (jKleine et al. 

2005 

Qin 

much stronger than observations. The key point 

et al. 2008 Kruijer et al. 

2014), and CAIs are 


is that, the dust grains also grow through coagula¬ 
tion, and as they grow, they settle to the disk mid¬ 
plane ( Dnllemond fc Dominik|2QQ5 ; Akimkin et ah 


2013), where they continue to grow bigger and 


get assembled into planetesimals that may even¬ 


tually grow to tens of kilometers in size (Stepin- 
ski fc Valageas|1997 ), possibly accompanied by in¬ 


ward drift ( Laibe et al.|2Q12 ). The majority of the 


oxygen and carbon originally frozen onto the dust 
grains in the atmosphere thus ends up inside solid 


bodies in the midplane of the disk ( 

Bergin et al. 

2010 Hogerheijde et al.||2011 

Najita et al.||2013). 


generally considered the oldest bodies of the so¬ 
lar system (Amelin et ah 2002). If systems such 


as TW Hya are similar to the solar system, then 
their planetesimals formed at an early stage will 
be rich in carbon- and oxygen-bearing ices. 

F.D. and E.A.B. are supported by grant 
NNX12A193G from the NASA Origin of Solar 
Systems Program. We thank Use Cleeves and 
Fred Ciesla for discussions. 
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